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Abstract 



In the bidirected minimum Manhattan network problem,, given a set T of n terminals in the 
plane, we need to construct a network N{T) of minimum total length with the property that the 
edges of N{T) are axis-parallel and oriented in a such a way that every ordered pair of terminals 
P^ , is connected in N{T) by a directed Manhattan path, bi this paper, we present a polynomial 

factor 2 approximation algorithm for the bidirected minimum Manhattan network problem. 
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1 Introduction 



A rectilinear network N = {V,E) in M? consists of a finite set V of points and horizontal and 
vertical segments connecting pairs of points of V. The length of A^ is the sum of lengths of its edges. 
'^ , Given a finite set T of terminals in M^, a Manhattan network [10] on T is a rectilinear network 

l/^ [ N(T) = {V, E) such that T C y and for every pair of points in T, the network N{T) contains a 

shortest rectilinear (i.e., Manhattan or li-) path between them. A minimum Manhattan network 

^^^ ^ on r is a Manhattan network of minimum possible length and the Minimum Manhattan Network 

C^ ' problem (MMN problem) is to find such a network (for an illustration, see Fig. 1). Note that there 

is always a minimum Manhattan network lying in the grid T(T) defined by the terminals (consisting 
of all line segments that result from intersecting each horizontal and vertical lines through a point 
inT). 

In this paper, we consider the following oriented version of the MMN problem. In the Bidirected 

5^ I Minimum Manhattan Network problem (which we abbreviate BDMMN problem), given a set T of 

terminals in the li-plane, we want to select a minimum-length subset N{T) of edges in the grid 
r(r) and to orient each edge in N(T) such that each ordered pair of terminals is connected by a 
directed Manhattan path (for an illustration, see Fig. 2). This oriented version of the minimum 
Manhattan network problem was formulated in [11] by M. Benkert and the second author of this 
note. Further we will assume that T does not contain two terminals on the same horizontal or 
vertical line, otherwise the problem does not have a solution. 

The minimum Manhattan network problem has been introduced by Gudmundsson, Levcopoulos, 
and Narasimhan [10]. Gudmundsson et al. [10] proposed an 0(n'^)-time 4-approximation algorithm, 
and an 0(n log n)-time 8-approximation algorithm. They also conjectured that there exists a 2- 
approximation algorithm and asked if MMN is NP-complete. Chin, Guo, and Sun [4] recently 
established that the problem is indeed NP-complete. Kato, Imai, and Asano [12] presented a 2- 
approximation algorithm, however, their correctness proof is incomplete (see [1]). Benkert, Wolff, 





Figure 1: A minimum Manliattan network Figure 2: A minimum bidirected Manhattan network 



Shirabe, and Widnrann [1] described an 0(n log n)-tinie 3-approxiination algorithm and presented 
a mixed-integer programming formulation of problem. Nouioua [16] and Fuchs and Schulze [7] 
presented two simple 0{n log ?i)-time 3-approximation algorithms. The first correct 2-approximation 
algorithm (solving the first open question from [10]) was presented by Chepoi, Nouioua, and Vaxes 
[3] and is based on a strip-staircase decomposition and a rounding method applied to the linear 
program from [1]. In his PhD thesis, Nouioua [16] described an 0(n log n)-time 2-approximation 
algorithm based on the primal-dual method. In 2008, Guo, Sun, and Zhu [8, 9] presented two 
combinatorial factor 2 approximation algorithms with complexity 0{in?) and O(nlogn) (see also 
the PhD thesis [19] of Schulze for yet another 0(n log n)-time 2-approximation algorithm). Seibert 
and Unger [18] announced a 1.5-approximation algorithm, however the conference format of their 
paper does not permit to understand the description of the algorithm and to check its claimed 
performance guarantee (a counterexample that an important intermediate step is incorrect was 
given in [7, 19]). Finally, a factor 2.5 approximation algorithm for MMN problem in normed planes 
with polygonal unit balls was proposed in [2]. 

Gudmundsson et al. [10] introduced the MMN problem in connection with geometric spanners. 
A geometric network A^ is a c-spanner (c > 1) for T if for each pair tj,tj € T, there exists a 
(tj,tj)-path in N of length at most c • ||tj — tj\\. In the Euclidean plane, the unique 1-spanner 
of T is the complete graph on T. In the rectilinear plane, the points are connected by several 
Manhattan paths, and the problem of finding the sparsest 1-spanner becomes non trivial. Minimum 
Manhattan networks are precisely the optimal 1-spanners. Analogously, the bidirected minimum 
Manhattan networks can be viewed as optimal 1-spanners of the directed grid F (T) obtained from 
F(r) by replacing each edge by two opposite directed arcs. Alternatively, bidirected Manhattan 
networks are roundtrip 1-spanners sensu [17] for the grid F (T) and for the set T of terminals. In 
both reformulations of bidirected Manhattan networks as directed 1-spanners of F (T) we suppose 
that, like in Manhattan Street Networks [5, 14, 20], an edge of F(T) participating in the resulting 
spanner can be directed only in one sense. Geometric spanners have applications in network and 
VLSI circuit design, distributed algorithms, and other areas [6, 15]. Lam, Alexandersson, and 
Pachter [13] used minimum Manhattan networks to design efficient search spaces for pair hidden 
Markov model alignment algorithms. 

In this paper, we propose a factor 2 approximation algorithm for the minimum bidirected Man- 
hattan network problems. We also solve in the negative Problem 6 from the collection [11] asking 



whether a specially constructed bidirected Manhattan network N0{T) is always optimal. 

Our algorithm employs the strip-staircase decomposition proposed in our previous paper [3] 
and subsequently used in other approximation algorithms for MMN. First we notice that each strip, 
oriented clockwise or counterclockwise, belongs to any bidirected Manhattan network. Then we show 
that all strips constituting larger sets, called blocks, have the same orientation. Since the strips 
from different blocks do not overlap, one can suppose that the algorithm orients the strips in the 
same way as in an optimal solution. Since the bases of a staircase and the strips touching it belong 
to a common block, it remains to construct in each staircase a completion of the strip's solution. 
Any optimal completion can be subdivided into two subnetworks which, forgetting the orientation, 
can be viewed as variants of the MMN problem for staircases. Such optimal (undirected) networks 
can be computed in polynomial time by dynamic programming. The algorithm then constructs a 
directed version of these networks by directing their edges and replacing some edges by two shifted 
oriented copies. 

We conclude this section with some notations. For a point p € M^ we will denote by p^ and 
py its two coordinates. For two points p, q of M^ we will denote by R{p, q) the smallest axis- 
parallel rectangle containing p and q (clearly, p and q are two opposite corners of R{p,q)). Let 
r = {ti, . . . ,t„} denote the set of n terminals (recall that T does not contain two terminals on the 
same vertical or horizontal line). For two terminals tj,tj E T we will set Rij = i?(ti,tj). We will 
say that the rectangle Rij is empty if Rij OT = {tj,tj}. Finaly, let F^ be the set of all ordered 
pairs {i,j) such that Rij is empty. 

2 The counterexample 



We start with an example showing that the bidirected network N0{T) defined in [11] is not optimal. 
N0(T) is defined in the following way: go through all empty rectangles Rij and orient the edges 
on the boundary of Rij clockwise if the line segment tjtj has positive slope and counterclockwise 
if tjtj has negative slope. A^0(T) is always a bidirected Manhattan network. In Fig. 3 we present 
an optimal bidirected Manhattan network (its length is 32) for a set of 5 terminals. For the same 
set of terminals, the length of the bidirected Manhattan network N0{T) presented in Fig. 4 is 34: 
N0{T) also includes the two dotted edges of the staircase not included in the optimal solution. 
Analogous larger examples show that the ratio between the length of N0 (T) and the optimum can 
be arbitrarily large. 
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Figure 3: An optimal bidirected Manhattan 
network 



Figure 4: The bidirected Manhattan network 

N0{T) 



3 Strips and staircases 

In this section, we briefly recall the notions of strips and staircases defined and studied in [3]; we 
refer to this paper for proofs and some missing details. An empty rectangle Rij is called a vertical 
strip if the x-coordinates of t j and tj are consecutive entries of the sorted list of all x-coordinates of 
the terminals. Analogously, a empty rectangle Rij is called a horizontal strip if the y-coordinates 
of tj and tj are consecutive entries of the sorted list of all y-coordinates of the terminals. The sides 
of a vertical (resp., horizontal) strip Rij are the vertical (resp., horizontal) sides of Rij- Notice 
that two points tj,tj may define both a horizontal and a vertical strip. We say that the rectangles 
Ri^i' and Rjj' form a crossing configuration if they intersect and they have the same slope. The 
importance of such configurations resides in the following property whose proof is straightforward: 

Lemma 3.1 If the rectangles Ri^i' and Rjji form a crossing configuration, then from the two couples 
of directed l\-paths connecting tj with tj/ and tj with tji one can derive two couples of directed l\- 
paths connecting tj with tj/ and tj with tj/. 
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Figure 5: A crossing configuration 



For a crossing configuration defined by the strips Ri^ii,Rjji we can define a rectangle with four 
tips as illustrated in Fig. 5. Let o and o' be the two opposite corners of this rectangle labeled in 
such a way that the four tips connect o with tj, tj and o' with tj/, tj/. Additionally, suppose without 
loss of generality, that tj and tj belong to Qi(o), i.e., to the first quadrant with respect to the origin 
o. Then tj/ and tj/ belong to Qsio'). Denote by Tjj the set of all terminals t^ G (T\{tj,tj})n Qi(o) 
such that (i) Ritk,o) DT = {tfc} and (fi) the region {q G Q2(o) : q^ < t^} U {q £ Q^io) : q"" < t%} 
does not contain any terminal of T. When Tjj ^ 0, we define the staircase 5j jij/ j/ as the union 
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of rectangles |J{i?(o', t^) : tk € Ti,j} \ -R(o, o'); see Fig. 6 for an illustration. The point o is called 
the origin of this staircase. Analogously one can define the set Tj/.j/ and the staircase 5j/ j/u ,■ with 



I, I, J belong to the second quadrant with 



origin a'. Two other types of staircases will be defined if tj,t 
respect to o and tj/, tj/ belong to Q4(o'). In order to simplify the presentation, further we will prove 
all results under the assumption that the staircase is located in the first quadrant. By symmetry, 
all these results also hold for the other types of staircases. 

Let a be the leftmost highest point of the staircase ^j jij/ j/ and let /3 be the rightmost lowest 
point of this staircase. By definition, Sijuiji DT = Tij. By the choice of Tjj, there are no terminals 
of T located in the regions {q G Q2(o) : q^ < a^} and {q £ Q^i^o) : q'-^ < f3^}. 





Figure 6: Staircase Sij^i^j 



Figure 7: Sij\i,j, with FfJ" and T°f oriented 




Figure 8: An example of a strip-staircase decomposition 



4 Blocks and generating sets 

For a strip Rij, the terminals tj and tj can be connected by two directed Manhattan paths only 
if we take the boundary of the rectangle Rij in the solution and orient it accordingly. Therefore 
the boundary of each strip Rij belongs to all bidirected Manhattan networks. Rij has only two 
orientations (clockwise and counterclockwise) producing the two directed Manhattan paths between 
tj and tj. We say that (the orientations of) two rectangles Rij and Ri'j' are compatible if they have 
the same slope and are oriented in the same way or if they have different slopes and are oriented 
in opposite ways. Clearly, two strips sharing an edge e of T{T) must be compatible. We extend this 
property to larger sets, called blocks. 

Let P C T be a maximal by inclusion set of terminals such that for all tj S P (the same) 
two opposite quadrants centered at tj are empty, i.e., their intersections with T consist only of tj; 
suppose that these empty quadrants are the second and the fourth quadrants Q2(ti) and Q4(ti). 
Now, let the points of P be sorted by x-coordinate. The ith block is the set of all terminals of T 
contained in the axis-parallel rectangle spanned by ith. and (i + l)th points of P, the first block 
consists of all terminals located in the third quadrant defined by the first point of P, and the last 
block consists of all terminals located in the first quadrant defined by the last point of P. From the 
definition follows that two terminals defining a strip belong to a common block (which some abuse 



of language, we will say that the strip itself belongs to this block). 

Lemma 4.1 In any bidirected Manhattan network N{T), all strips from the same block B are 
compatible. 

Proof. Consider a graph whose vertices are the strips from B and two strips are adjacent if and 
only if their boundaries share an edge of T{T) or a terminal. This graph is connected because its 
two subgraphs induced by horizontal and vertical strips are connected and any terminal of B defines 
in 5 a vertical and a horizontal strip which are adjacent in this graph. Therefore, if B contains 
incompatible strips, then we can find in this graph two adjacent incompatible strips Rij and Rj'^k- 
Since two strips sharing an edge are compatible, Rij and Rj'^k necessarily share a terminal, say 
j = j' . We can suppose without loss of generality that tj € Qi{tj) and t^ G Qsi^j), otherwise the 
boundaries of these strips overlap. Hence tj ^ P, and Q2(tj) or (54(tj) is not empty. Suppose that 
there is a point t; € Q4{tj) (see Fig. 4 for an illustration). Since Rij and Rj^k are incompatible, 
tj is the head or the tail of both edges incident to tj in Q4{tj), say the tail. Since t; G Q^itj), the 
rectangle Rj^i is also included in Q4(tj). Since any (directed or not) Manhattan path between tj 
and t; is included in Rji and both edges of N{T) D Rj^i incident to tj are directed away from tj, we 
will not be able to produce a directed Manhattan path from t; to tj in N{T), a contradiction. D 



-® 



Figure 9: To the proof of Lemma 4.1. 

We continue by adapting to BDMMN the notion of a generating set introduced in [12] for MMN 
problem: a generating set is a subset F of ordered pairs of terminals of T with the property that 
a directed subnetwork of T{T) containing directed Manhattan paths between all pairs of F is a 
bidirected Manhattan network for T. Let F' be the set of all ordered pairs of terminals defining 
strips. Let also F" be the set of all ordered pairs (tj/,t/) and (t/,tj') such that there exists a 
staircase <Siju'j' with t; belonging to the set Tij of all terminals defining the corners of Sijuiji. 
The proof of the following result closely follows the proof of Lemma 3.3 of [3]. 

Lemma 4.2 F := F' U F" is a generating set. 

Proof. The set F0 of empty pairs is clearly a generating set. Let A^ be a bidirected rectilinear 
network containing directed ^i-paths for all ordered pairs in F. To prove that A'^ is a bidirected 
Manhattan network on T, it suffices to establish that for any arbitrary pairs {k, k'), {k' , k) G F0 \F, 
in N there exists a directed Manhattan path from t^ to tk' and a directed Manhattan path from 
tfc/ to tjt- Assume without loss of generality that t^, < t^ and t|, < t|. Since {k,k') G F0, the 
rectangle Rk,k' is empty. The vertical and horizontal lines through the points t^ and tk' partition the 
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Figure 10: To the proof of Lemma 4.2 



plane into the rectangle Rk,k'^ four open quadrants and four closed unbounded half-bands labeled 
counterclockwise 'Hi,'H2,'H3, and Ti^ (see Fig. 10). Consider the leftmost terminal tj^ of Hi (this 
terminal exists because t^ € "Hi)- Now, consider the rightmost terminal tj/ of "Ha such that t^ < tf 
(again this terminal exists because t^/ € Tis and t^, < if ). By the choice of tj^ and tj/ , the rectangle 
i?jj j' is the leftmost vertical strip crossing the rectangle Rk,k'- Analogously, by letting tj2,tj^, and 
tj2 be the rightmost terminal of Ti^, the lowest terminal of 7^4 and the highest terminal of T-L2, 
respectively, we obtain the rightmost vertical strip Ri^^i' , the lowest horizontal strip Rj-^j' , and the 
highest horizontal strip Rj^j; crossing the rectangle Rk,k'- Notice that the strips Rj^.j' and Ri2,i' 
as well as the strips Rjj^j' and -Rj^^j' constitute crossing configurations. 

Now, we will prove that N contains a directed ^i-path from tjj to t^ and a directed ^i-path from 
tfc/ and tj-^ . We distinguish three cases. If t^ = t j/ , then Ri^^k = Ri2,i' is a strip and thus (fc, 12) & F. 
If tk = tji^, then the strips Rj2,k and Ri^^'^ form a crossing configuration. By Lemma 3.1, from the 
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directed ^i-paths of A^ running from tj^ to t^ and from tjj to tj/ , we can derive a couple of directed 
Zi-paths connecting t^ with tjj. Finally, if t^ ^ {tj/ ,t,/}, we assert that the crossing configuration 
Ri2,i' and Rj2,j' defines a staircase Si' j'\i2,J2 such that t^ belongs to T^i ji . Indeed, let o be the 
highest leftmost intersection point of the strips Ri2,i' and Rj2,j' (see Fig. 10). Since i?(tfc,o) is 
contained in the empty rectangle Rk,k', we conclude that i?(tfc,o) n T = {t^}. Moreover, by the 
choice of t jj and tj^ , the unbounded half-bands {q G T-L-^ : q^ > t^i } and {q € 7^2 : <f > ^j' } do 
not contain terminals (in Fig. 10, the shaded region does not contain terminals), thus establishing 
our assertion. This implies that tk € Tj^i ji , whence {k,i2), {12, k) € F. Therefore, in all three cases 
the terminals t^ and tjj are connected in A^ by a couple of directed /i-paths. Using a similar 
analysis, one can show that t^/ and tj^ are also connected in A^ by a couple of directed ^i-paths. By 
construction, the rectangles Rk,i2 and Rk'ji form a crossing configuration and thus, by Lemma 3.1, 
there is a couple of /i-paths of A^ between the terminals t^ and t^/, concluding the proof. D 



For a staircase 5j j|j/ j/, 



let T^- be the set consisting of Tij , the four terminals t j , t j , t j/ , tj/ of the 



bases of Sij^iijt, and the terminals defining strips touching the boundary of Sij^itji. 

Lemma 4.3 T- ■ is included in a block. 

Proof. Suppose by way of contradiction that there exists t^ G P such that two terminals t; and 
tm of Tj"t- belong to different quadrants Qi{tk) and Qsl^fc) (and the quadrants (52(tfc) and Q4{tk) 
are empty). Since the interiors of Siju/j/, of the rectangles Riy,Rjj', and of the regions Q' , Q" 
defined in previous section are all empty, t^ can be located only in the first quadrant defined by a 
concave vertex of Siju'j'- But in this case, we can find two terminals of Sijuji located in Q2{^k) 
and Q4(tfc), contrary to the choice of t^ in P. D 



5 The algorithm 

By Lemma 4.1, all strips of any block are compatible. Since the strips from different blocks are 
edge-disjoint, the algorithm can test the two possible orientations of each block independently of 
the orientations of other blocks. Thus, we can suppose that the algorithm oriented the strips of each 
block in the same way as in an optimal bidirected Manhattan network N*(T). Let N'(T) be the 
union of all boundaries of strips directed this way. By Lemma 4.3, the bases of a staircase Sijuiji 
and the strips touching Siju'ji belong to the same block B, therefore they are all compatible and 
their orientation can be assumed fixed. Since the bases of Siju/j/ have the same slope, they are 
oriented both clockwise or both counterclockwise. From [3] we know that any strip may touch the 
boundary of a staircase but cannot cross its interior. Let N- ■ be the intersection of Sj_ju/j/ with 
the union of the boundaries of the strips from B, i.e., N^ • = N'{T) D Siju/ji. Hence, by Lemma 4.2 
it remains, for each staircase Siju/j', to complete A^| ■ to a local bidirected network N-'^ such that 
any pair (tj',t/), (t/,tj/) with t; G Tjj can be connected in N", U N'{T) by a directed Manhattan 
path. This must be done in such a way that the length of the network Nij = N"- \ N^ ■ is as small 
as possible. Let N* ■ be an optimal completion of N[-. 

Suppose that i?j j/ and Rjj' are the vertical and the horizontal bases of 5j ju/ j/ (see Fig. 6). 
Let R'^^i = Ri^i' n 5jj|j/j/ and R'- ■, = Rjj' nSij^i/ji. Suppose that in algorithm's and optimal 
solutions, these strips are oriented clockwise. Hence the leftmost vertical side Sj/ of R'^ ^, is upward, 
the opposite side Si is downward, the upper horizontal side Sj of R'- ■, is to the right, and its opposite 
side Sj' is to the left. For a terminal t € Tij, denote by t and t the tail and the head of the 
directed edges of N-j incident to t. Set Tij = {t : t G Tij} and Tij = {t : t e Tij} (they are 
all vertices of the grid T{T)). 

Any optimal completion iV* • of N^ ■ can be decomposed into two edge-disjoint subnetworks N*ij 
and N*ij, such that A^ ij contains the edges on the directed Manhattan paths running from Sj' Usj 
to the points of Tij and N*ij contains the edges on the directed Manhattan paths running from 
the points of lij to Si U Sj/. The length of N*ij cannot be smaller than the length of an optimal 
(non-oriented) network Aij in T(T) connecting the points of T ij to Sj' U Sj by Manhattan paths. 
Analogously, the length of N*ij cannot be smaller than the length of an optimal fnon-oriented) 
network Bij in T{T) connecting the points of T ij to Sj U Sji. At the difference of N*ij and N*ij, 
the networks Aij and Bij are not edge-disjoint. However, we can compute optimal Aij and Bij 
in polynomial time using dynamic programming because each of the sets of terminals Tij and 



Tij also constitute a staircase and to compute an optimal solution we will have to solve only a 
polynomial number of subproblems (this problem is similar to Steiner arborescence or to MMN 
problems for terminals on a staircase). 

The algorithm computes by dynamic programming an optimal network Aij for connecting T ij 
to Sj' U Sj and an optimal network Bij for connecting Tjj to Sj U Sj'. The dynamic programming 
constructs Aij recursively in the following way: there exists a point t € Tij which is either 
connected in Aij to s^/ by a horizontal segment s' or to Sj by a vertical segment s", say the first. 
Then the problem is subdivided into two smaller subproblems, one for the points of T ij located 
strictly above s' and the union Sj' U s' and another for Sj/ U Sj and the points of T ij located strictly 
below s' (the case when t is connected vertically is analogous). The construction of Bij is similar 
(see the first two networks in Fig. 5 for an illustration of Aij and Bij). 

Finally, the algorithm "rounds" each of the networks Aij and Bij in order to produce directed 
networks Aij and Bij having lengths at most twice the lengths of Aij and Bij, respectively (see 
the last two networks from Fig. 5). The algorithm returns Nij = Aij U Bij as a local completion 
of N--. The networks Aij and Bij are constructed in the following way. Let vq,vi, . . . ,Vk-i 
and hi,h2, ■ ■ ■ ,hi be the consecutive horizontal and vertical lines of the grid T{T) intersecting 
the staircase Sij^i/j/ and numbered in such a way that Si' C vo,Si C vi and Sj' C hi,Sj C /i2. 
(If the bases are oriented counterclockwise, then we consider the same lines but we index them 
vi,V2, ■ ■ ■ ,Vk and ho, hi, ... , /i/_i.) Let Ffj^" (respectively, T"j^) be the subgrid of r(r) n Sij^i/jr 
induced by vertical and horizontal lines with even indices (respectively, with odd indices). (Notice 
that if we orient the horizontal edges of Tf^^"^ to the right and the vertical edges upward, and the 
horizontal edges of ^"'¥' to the left and the vertical edges downward, then we obtain a network 
which is well-known in the literature as a Manhattan Street Network (see Fig. 7) [5, 14, 20].) 

The algorithm recursively derives Aij from A^j and Bij from Bjj. We show how to construct 
Aij, but each step of the algorithm is performed for both Aij and Bij (before going to the next 
step). First, in Step 1, we insert in Aij the edges of T^""^^ which have their support in Aij and 
orient them upward or to the right. The remaining directed edges are added in order in which the 
segments of Aij have been added by the dynamic programming algorithm. For a current t ^ let s 
be the vertical or horizontal segment of Aij connecting t to the previously defined part of A ij . If 
s belongs to a horizontal or vertical line with odd index m, say s belongs to h^, then in Step 2 we 
include in A^j the segments s', s" oriented to the right which correspond to paths of Tf^^"' obtained 
by intersecting /im-i and hm+i with the vertical lines passing via the ends of s. If Sj C /im-i (i-e., 
?n = 3), then s' is not added. Additionally, we remove from A^j all vertical edges whose lowest 
end-vertex is comprised between s' and s" (by the construction, this pruning operation removes the 
edges that are no longer used by directed Manhattan paths in Aij). In Step 3, we proceed the 
points of lij n T°'^ in the same order as in Step 2. Let c be the vertical segment connecting t to 
hm-i- If c does not belong to Bij after Step 2 (Fig. 13, Step 3 (a)), we add c oriented upward to 
Aij. Otherwise, we consider the cell of r(r) whose boundary contains c and a subsegment b' of s' 
(Fig. 13, Step 3 (b)), we remove b' and add to Aij the alternative path around this cell consisting 
of an upward twin c' of c and a twin h of b' oriented to the right (6 is a subsegment of s). We call 
such a path a replacement path. Let N'^'- := N[ ■ U Aij U Bij and let 1{N['-) be its length. 

We conclude this section with the pseudocode of the algorithm constructing the directed networks 
Aij and Bij from Aij and Bij: 



®— ' 



n 



SH-I 



13 



1 1 » ^ 



(a) Ai,, 



(b) B.,, 



r^ 



iL " 



S^ 



> D> » 




> □ » <^ 






QVKS) 



SH-I 



1 < ! » » 



I 



(c) ~ii. 



(d) 



Figure 11: The networks Aij,Bij and the rounded directed networks Aij, B 
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Algorithm 1 Construction of Aij U i?j^ 



Require: Networks Aij and -Bjj-- 

Step 1: Insert in Aij the edges of T?'!'^'^ which have their support in Aij and orient them upward 
or to the right. 

Step 2: For each remaining directed edge (in order in which the segments of Aij have been 
added by the dynamic programming algorithm): 

If Sj (^ hrn^-i (i.e., 771 7^ 3), then include in Aij the segment s' oriented to the right. 

Include in Aij the segment s" oriented to the right. 
Remove from Aij all vertical edges whose lowest end- vertex is comprised between s' and s". 
Perform the Steps 1-2 for computing Bij from Bij. 
Step 3: For each point of T ij n F"'^'^ (proceeded in the same order as in Step 2): 

If c does not belong to Bij, then add to Aij the segment c oriented upward. 

Otherwise, remove b' and add to Aij the alternative path consisting of c' oriented upward 
and b oriented to the right (see Fig. 13). 
Perform the Step 3 for each point of j j ,■ n Tf^f"'. 

6 The analysis of the algorithm 

Now, we will show that the algorithm described in previous section returns a bidirected Manhattan 
network and that the length of this network is at most twice the length of an optimal bidirected 
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Figure 13: Steps 2 and 3 of Algorithm 1 



Manhattan network. 

Lemma 6.1 The supports of Aij and Bij are disjoint. Moreover, in Aij there exists a directed 
Manhattan path from Si' U Sj to each point of T ij and in Bij there exists a directed Manhattan 
path from each point of lij to Si U Sj'. 

Proof. By the algorithm, after Step 2 the supports of Aij and Bij are disjoint, however these 
networks are not yet feasible. Step 3 ensures feasibility of Aij (and Bij) by connecting each 
terminal t E T ij to the network A^j computed in Step 2, using either a vertical segment c (Fig. 
13 Step 3 (a)) or a replacement path {b,c'} (Fig. 13 Step 3 (b)). We will prove now that after 
Step 3 the networks A^j and Bij are disjoint and feasible. 

By construction, Aij connects each terminal of T ij to Si' Usj by a Manhattan path. Therefore, 
it suffices to show that this property is preserved each time when we orient a new part of the 
network Aij. This is obviously true when we orient all edges of F?';'-^" which have their support 
in Aij. Now, suppose that a segment s added in Aij to connect a terminal t is replaced by two 
directed segments s' and s". Then all vertices of Tij connected in A^j via s (they are all located 
above s) will be now connected by directed Manhattan paths going via s". On the other hand, t 
is connected in Aij via c and s' if c does not belong to Bij and via the replacement path {b,c'} 
and s' otherwise. We assert that in the last case, c' and b will not be used by Bij. Indeed, due 
to the pruning operation in Step 2, each t is incident to at most one outgoing edge of Bij. Now, 



■(— 



since we used a replacement path, the segment c of T°j is still in Bij after the pruning operation 
in Step 2. This shows that b cannot belong to Bij. On the other hand, the end-points of c' do not 
belong to T ij, otherwise two distinct points of T ij U lij will lie on the same vertical or horizontal 
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line, which is impossible. Now, since c' belongs to Tf""^^, by the algorithm, c' can be involved in a 
replacement path of Bij only if it is incident to a point of J ij- Therefore c' does not belong to 
Bij either. D 

Lemma 6.2 l{Nij) = l(Aij U Bij) < 2 l{N*j) + l{Nlj), where N*j is an optimal completion of 

N' ■. 

Proof. By construction, the length of AijUBij is at most 2 l{Aij)-'r2 l{Bij) plus the total length 
A of the edges c or c' orthogonal to s' defined in the algorithm. Since any horizontal or vertical line 
crosses at most one such edge, A is at most 1{N[-), hence l{Nij) < 2 {l{Aij) + l{Bij)) + l[N[-). 

By optimality of Aij and S^j, l{Ai^j) + l{Bi^j) < /(iV^ij) + /(^ij) = KN*j). □ 

Let N{T) be the network obtained as the union of N'(T) and all local completions Nij = 
Aij U Bij taken over all staircases. Let N*{T) be an optimal solution of the BDMMN problem 
having N'(T) as a subnetwork. Then N*(T) can be viewed as the disjoint union of N'(T) with the 
local completions N*j = {N*{T) \ N'{T)) nSijii>ji of iV^. It was shown in [3] that the interiors of 
two staircases are disjoint. Since in the BDMMN problem there are no two terminals on the same 
horizontal or vertical line, two staircases cannot intersect in an edge, thus the intersection of two 
staircases is empty or a subset of terminals. Hence the local completions N*- are pairwise disjoint (as 
well as the local completions Nij). By their definition, the networks N- , are also pairwise disjoint. 
Using this disjointness property, Lemma 6.2, and summing up over all staircases, we obtain that 

1{N{T)) = l{N'{T))+Y,KN,,)<KN'{T))+Y,{2l{N*^) + liNl^)) 
< 2 liN\T)) + 2 1{N*{T) \ N'{T)) = 2 1{N*{T)), 

The time complexity of the algorithm for the BDMMN problem is dominated by the execution 
of the dynamic programming algorithm that computes Aij and Bij for every staircase Sijuij/- A 
staircase SijLrjr contributes 0{\Tij\^) to the total complexity of the algorithm. Since each terminal 
belongs to at most two staircases, the total complexity of the algorithm for the BDMMN problem 
is 0{n^), establishing the main result of the paper: 

Theorem 6.3 The network N(T) computed by the algorithm in 0{n^) time is a factor 2 approxi- 
mation for the BDMMN problem. 

Open question: Is the BDMMN problem polynomial or NP-hard? BDMMN reduces only to 
staircases, avoiding thus the difficulty occurring in the MMN problem due to the interaction between 
strips and staircases. However, in the directed version, for each staircase we have to compute two 
disjoint but not necessarily optimal directed networks Nij and Nij which together provide an 
optimal completion of N- ■ . 
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